A systematic approach to bicontinuous cubic phases in ternary amphiphilic systems 



U. S. SchwarzB and G. Gompper 
Max- Planck- Institut fur Kolloid- und Grenzflachenforschung, 
Kantstr. 55, 14513 Teltow, Germany. 
(February 1, 2008) 

' The Fourier approach and theories of space groups and color symmetries are used to systemat- 

, ically generate and compare bicontinuous cubic structures in the framework of a Ginzburg-Landau 

■ model for ternary amphiphilic systems. Both single and double structures are investigated; they 

correspond to systems with one or two monolayers in a unit cell, respectively. We show how and why 
single structures can be made to approach triply periodic minimal surfaces very closely, and give 
improved nodal approximations for G, D, FWP and P surfaces. We demonstrate that the relative 
stability of the single structures can be calculated from the geometrical properties of their interfaces 
only. The single gyroid G turns out to be the most stable bicontinuous cubic phase since it has the 
smallest porosity. The representations are used to calculate distributions of the Gaussian curvature 
CN \ and 2 #-NMR bandshapes for C(P), C(D), S, C(Y) and F-RD surfaces. 
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Amphiphiles are molecules which have both hydrophilic and hydrophobic parts. In a ternary mixture of amphiphiles, 
water and oil, their ability to stabilize water-oil interfaces leads to structure formation on the nanometer scale 
Depending on concentration and temperature, many different phases are found to be stable. One intriguing aspect 
£j ■ of amphiphilic polymorphism is the existence of ordered bicontinuous phases which can be traversed in any direction 
in both the water and the oil regions. Here macroscopic order has been demonstrated by the appearance of Bragg- 
peaks in diffraction patterns [Q] and bicontinuity by measuring diffusion properties with nuclear magnetic resonance 
|3|. Clearly any bicontinuous ordered phase has to have a three-dimensional Bravais-lattice. In experiments, cubic 
symmetry is observed in most cases The amphiphilic monolayers of bicontinuous cubic structures, which separate 

regions with water from regions with oil, are often modelled by triply periodic surfaces of constant mean curvature 
P,fl0[ . By changing experimental parameters like temperature or salt concentration, their preference to bend towards 
the water or the oil regions can be made to vanish. Then one can describe them by triply periodic minimal surfaces 
(TPMS). Similarly, the bilayers of bicontinuous cubic phases in lipid- water mixtures by symmetry have no spontaneous 
curvature, so that their mid-surfaces can also be modelled by TPMS Bicontinuous cubic phases have gained 

renewed interest recently, for example as space partitioners in biological systems |l2|| , as amphiphilic templates for 
mesoporous systems jH| and for the crystallization of membrane proteins p4fl . 

In this paper, we report on a systematic theoretical investigation of bicontinuous cubic phases in ternary amphiphilic 
systems with vanishing spontaneous curvature. We restrict ourselves to bicontinuous ordered phases with cubic 
symmetry since these are the ones which are mainly observed in experiments. It has been shown recently by Gozdz 
and Holyst |l5| ] in the framework of real-space minimization that various cubic bicontinuous phases can be generated 
as local minima of the Ginzburg-Landau model for ternary amphiphilic systems introduced previously by Gompper 
and Schick . The position of the interfaces between oil and water regions then follow as iso-surfaces to a scalar 
field $(r). Our analysis is based on the Fourier ansatz in combination with the theories of space groups and color 
symmetries; it allows to generate the relevant structures in a very systematic way and to reduce the representations 
to a relatively small number of variables. This in turn allows very efficient numerics and makes it easy to document 
and reuse the results. 

In order to analyze our results, we use an effective interface Hamiltonian, which is derived from the Ginzburg-Landau 
theory J!?]]. Due to the oil- water symmetry of the Ginzburg-Landau model, the spontaneous curvature vanishes. Since 
the elastic energy of interfaces is governed by their bending rigidity, the interfaces of the structures which have only one 
amphiphilic monolayer are very close to cubic TPMS. A similar approach, which employs a Ginzburg-Landau model 
for binary fluid mixtures, has been used in Refs. ]l~8| , |il| . However, since the elastic energy of interfaces is governed 
by their tension in this case, the numerical minimization is much less stable. In fact, our tabulation of the first few 
Fourier amplitudes for the four best-known structures G, D, I-WP and P provides a considerable improvement over 
the widely used nodal approximations for TPMS [^(j . Since most of the detailed analysis of the properties of cubic 
TPMS has been focused on these structures so far, we want to use our representations to characterize the structures 
C(P), C(D), S, C(Y) and F-RD in more detail. 
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Our analysis in terms of interfacial properties also allows us to investigate the relative stability of the various 
bicontinuous cubic phases in terms of their geometrical properties. Although the identification of a certain bicontinuous 
cubic structure is a tedious procedure best accomplished by combining several experimental techniques like neutron 
or X-ray small angle scattering, transmission electron microscopy and nuclear magnetic resonance, in most cases 
one has found the stable phases to be either gyroid or diamond structures. Here we derive an universal geometrical 
criterion for their relative stability in the case of vanishing spontaneous curvature of the interfaces which explains 
these experimental findings and is in excellent agreement with our numerical results. 

The paper is organized as follows. We begin by introducing the relevant concepts from differential geometry and 
crystallography and their applications to amphiphilic systems in Sec. [H| In Sec. [II we explain how to construct 
Fourier series and nodal approximations for amphiphilic structures. In Sec. [fy| we systematically generate, evaluate 
and compare various single and double structures by using the Ginzburg-Landau theory, as well as the interface 
Hamiltonian which follows from it. In Sec. [v| we use the representations obtained to comment on the experimental 
identification of bicontinuous cubic phases in amphiphilic and mesoporous systems. 



II. GEOMETRY OF BICONTINUOUS PHASES 



If we characterize a bicontinuous phase by its interfaces, the simplest possible case is that they form one triply 
periodic surface (TPS), i.e. a surface with a three-dimensional Bravais-lattice. Any TPS divides space into two 
unconnected but intertwined labyrinths. In ternary amphiphilic systems, the surface is covered with an amphiphilic 
monolayer and the two labyrinths are filled with water and oil, respectively. Therefore any structure corresponding 
to a TPS is bicontinuous in water and oil. The TPS is called balanced if there exists an Euclidean transformation a 
which maps one labyrinth onto the other. Oil and water have equal volume fractions in this case. 

In the mathematical context, triply periodic minimal surfaces (TPMS) have been investigated extensively |2l|j22"f . 
In order to locally minimize their surface area, minimal surfaces have vanishing mean curvature H = (I/-R1 + l/i?2)/2 
at every point on the surface, where R\ and R2 are the two principal radii of curvature. Therefore, their Gaussian 
curvature K = l/RiR 2 = —1/R\ is non-positive and the surface is saddle-shaped everywhere. Up to 1970 only five 
TPMS had been known by the 19th-century work of Schwarz and his students, with P, D and C(P) being cubic. Then 
Schoen described 12 more [p3[ , including the cubic cases G, F-RD, I-WP, 0,C-TO and C(D)- Their existence was 
proven by Karcher in 1989 P4|. More TPMS have been found by Karcher and Polthier j2j|, Fischer and Koch |2(| 
and others, but all of them seem to be more complicated than the ones known to Schoen. For cubic TPMS, exact 
representations are known only for P, D, G and I-WP. These Weierstrass representations originate from complex 
analysis and are explained in Appendix |A|. With their help the surface areas of these TPMS can be calculated exactly 
in terms of elliptic functions; the resulting values are given in Tab. |. 

Being balanced and being minimal are two independent possible properties of a TPS. From the cubic TPMS known 
to Schoen, F-RD, I-WP and 0,C-TO are not balanced and divide space into two labyrinths with unequal volume 
fractions. On the other hand, the TPS defined by cos a; + cosy + cosz = is balanced with a being a translation 
along half the body diagonal of the conventional unit cell, but far from minimal. If a TPS contains straight lines, it 
is called spanning |27| . It was shown by Schwarz that any spanning TPMS is symmetric with respect to a reflection 
around such a line. Since this operation exchanges the two labyrinths, any spanning TPMS is balanced. The only 
known cubic TPMS, which is balanced but not spanning, is the gyroid G, where the symmetry operation exchanging 
the two labyrinths is the inversion. 

For our purposes it is best to investigate not the local concentrations of oil and water separately, but their difference 
$(r) = po(j) — Pw{y)- Since we assume incompressibility, po( r ) + Pw{y) = 1 and the local concentrations of oil and 
water can be reconstructed from <f>(r). The amphiphilic monolayers now correspond to the $ = iso-surfaces. The 
interchange of oil and water amounts to $ — > — $. For balanced surfaces this is a symmetry operation when combined 
with the Euclidean operation a, that is $(ar) = — $(r). However, it is important to note that this is not a space 
group symmetry — which would only allow for a rigid motion, but not for a change of sign. Additional symmetries 
which arise when an Euclidean operation is combined with a permutation are called color symmetries; with only water 
and oil present we have the simplest case of a black and white symmetry [^8] [30|. When water and oil are colored 
white and black, respectively, the $ = iso-surface obtains differently colored sides and thus can be considered to be 
oriented. We assume the uncolored structure (i.e. the unoricnted surface) to have space group Q. Then all Euclidean 
operations which do not interchange the two labyrinths form a subgroup TL of Q of index 2. By definition TL is the 
space group of the colored structure (i.e. the oriented surface). The quotient group is isomorphic to the cyclic group 
of order 2, that is Q /TL = Z2 = {1, a}. In general, any balanced TPS is characterized by a group-subgroup pair Q — TL 
of index 2. For example, for the balanced TPS defined by cos a; + cosy + cosz = 0, we have TL = Pm3m (No. 221) 
and Q = Im3m (No. 229). 
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In order to describe a balanced TPS which does not intersect itself, certain restrictions apply for the extension from 
Ti to Q (for example, if a is a rotation, it has to be 2-fold). This has been used by Fischer and Koch to develop a 
crystallographic classification of all cubic TPMS |27|]2l| . In particular they showed that only 34 cubic group-subgroup 
pairs Q — Ti with index 2 can correspond to cubic balanced TPMS. Although there is no way to systematically list 
all balanced TPMS belonging to a possible group-subgroup pairing Q — Ti, Fischer and Koch argued that one can 
completely enumerate all cubic spanning TPMS which have straight lines forming a three-dimensional network whose 
polygons are spanned disc-like by the surface. They find P, C(P), D, C(D), S and C(Y), with the latter two described 
by them for the first time. In our work we systematically generate these structures together with the only known 
balanced but non-spanning cubic TPMS G and the non-balanced ones I-WP and F-RD. Examples for spanning TPMS 
are shown in Fig. |[ 

One also can consider structures with two TPS. For example, parallel surfaces can be constructed on either side of 
any given TPS. Both parallel surfaces are TPS, each of which has two labyrinths, which are topologically equivalent to 
those of the initial structure. However, with respect to the overall structure, these labyrinths are now filled with the 
same component and separated by a bilaycr which is filled with the other component. One might call this structure 
tricontinuous, but we rather prefer to use the term double structures for structures of two TPS — in contrast to single 
structures containing one TPS. 

Experimentally, little is known about the exact structure of bicontinuous cubic phases in ternary amphiphilic 
systems. One exception is the system DDAB-water-styrene, for which five different double structures have been found 
to be stable, each of which can be regarded as an oil-filled bilayer draped onto a different TPS (3l). However, with 
microemulsion phases being so ubiquitous for these systems, single structures can be expected also to be stable, 
and in fact have been reported in Ref. |32| ]. In contrast to the situation for (almost balanced) ternary amphiphilic 
systems, the structure of bicontinuous cubic phases is well-established for many binary amphiphilic systems [^]-^|. 
Here, mostly inverse phases are found, in which the two labyrinths are filled with water and the TPS is covered with 
an amphiphilic bilayer. These structures can be considered to be double structures of a ternary system where the oil 
has been removed from the bilayer. However, if the thickness of the bilayer is small compared to the lattice constant, 
it seems more appropriate to model them by single structures of a ternary system, where the water regions on both 
sides of the bilayer are distinguished artificially by labelling them 'water I' and 'water II', respectively (B3|. Thus our 
results for single structures can also be applied to bicontinuous cubic phases in amphiphile- water mixtures. 



III. FOURIER ANSATZ 



In order to construct a Fourier ansatz for a given group-subgroup pair Q — Ti, we need a Fourier ansatz for space 
group Ti, together with the implementation of a certain black and white symmetry. In ternary amphiphilic systems, 
this corresponds to a balanced single structure. For non-balanced single structures and double structures there is no 
color symmetry and one only needs the Fourier ansatz for a given space group Ti = Q . 



It is well known how to choose the Fourier ansatz for a given space group Ti [p4 35 ■ In the cubic system the 
translational properties of Ti correspond to one of the three cubic Bravais lattices simple cubic (P), body centered 
cubic (I) or face centered cubic (F). The Fourier series then reads 

$(r) = ]T 4>k e lK r (1) 

{K} 

where {K} is the reciprocal lattice and the are the (complex) Fourier amplitudes. However, not all of them are 
independent. Since the density field $(r) has to be real, it follows from Eq. (|l|) that $k = ^-k- With $k = A^ — iBm 
, Eq. ([!]) becomes 

$(r) = i A * cos(K • r) + B K sin(K • r)} . (2) 

{K} 

Here and Bk are real and can be obtained from $(r) as 

A K = - [ dv $(r) cos(K • r), B K = - [ dr $(r) sin(K ■ r), (3) 
v J v J 

where the integration runs over some unit cell of volume v. Thus, for even or odd functions, -Bk = or = 0, 
respectively. 

Further restrictions on the Fourier amplitudes arise when we consider symmetry operations r — > Pr + 1, which 
are not pure translations. Here, P denotes a rotation, inversion or roto-invcrsion, and t a translation which does not 
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belong to the Bravais lattice. Thus, t is non- vanishing for screw axes and glide planes. It now follows from Eq. (|lj) that 
(f>p-i K = <I>k e lK t . Let H be the point group of the given space group TL. To each operation in H there corresponds 
a pair (P, t) in TL. We group all reciprocal lattice vectors, which are related by an operation of H, into so-called 
Fourier stars. All members of a given star have the same wavelength. Only few stars will have the same wavelength 
and in most cases there is only one star with a certain wavelength. Therefore, the stars can be ordered according 
to decreasing wavelength and numbered with the index i. For each star we arbitrarily choose one representative Kj. 
With $k, = Ai — iBi, Eq. (|l|) becomes a sum over Fourier stars, 

i 

where the geometric structure factors are defined by 

A Ki (r) = cosKP^KO • r + K, ■ t], (5) 

H 

S K< (r) = MiP^Ki) ■ r + Ki ■ t] . (6) 

H 

Their precise form is obtained by explicitly enumerating all non-translational symmetry operations of space group 
TL. The results of this procedure can be taken from Ref. |}4j if they are corrected by a factor 1/rriL, where mj, is 
the multiplicity of the Bravais lattice (1, 2 and 4 for P, I and F, respectively). In Eq. (S), n denotes the order of H 
and rii the multiplicity of Kj. If H does not contain the inversion, we can extend the definition of a Fourier star by 
using $k = $-k- This does not change Eq. (0) except that now we have to replace H by the so-called Laue group, 
that is H extended by the inversion. All structures investigated in this work have m3m with order n = 48 either as 
their point group or as their Laue group, thus the members of one star are generated by all possible permutations 
and changes of sign of the components of Kj. With A4 = A^. and B± = B^., every Fourier star has two variables. 
This reduces to one when $ is an even or odd function of r. 

We now address the question how to implement the additional black and white symmetry. Since TL has index 2 in 
Q, it follows from the theorem of Hermann |29) that TL has either the same point group or the same Bravais lattice as 
Q. If TL and Q have the same point group, their Bravais lattices have to be different. Thus the Euclidean operation a, 
which maps oil regions onto water regions, has to be a translation t a which transforms one cubic Bravais-lattice into 
another. There are only two possibilities: for t a = a(x + y + z)/2 a P-lattice becomes a I-lattice, and for t a = ax/2 
a F-lattice becomes a P-lattice. The operation $(r + t a ) = — $(r) in combination with Eq. ([IJ) leads to additional 
reflection conditions: h + k + l = 2n+l for P — ► I, and h,k,l = 2n + l for F — > P. Note that these reflection conditions 
are complementary to the well-known ones for I and F (h + k + l = 2n and h + k,h + l,k + l = 2n, respectively). Thus 
in the case of identical point groups, the form of the Fourier stars in Eq. ([|) for space group TL stays the same, but 
the sum now runs over a reduced set of Fourier stars. 

If TL and Q have the same Bravais lattice, their point groups have to be different. Thus the Euclidean operation 
a has to be a point group operation P a which extends one cubic point group into another one. There are five cubic 
point groups and six ways to extend one of them into another by some P a . From $(P Q ,r) = — $(r) and Eq. (|l]) we 
now find & P -i K = — < I > Kexp(iK • t Q ). Thus, different Fourier stars merge to form new ones. However, in all cases 
considered in this work there is no need to consider new Fourier stars, since we only deal with the simple case that 
P a is the inversion. Then t Q = and one finds Ai = for all i. The same result of course follows from Eq. (Q), since 
now $ is an odd function of r. 

Tab. H lists the group-subgroup pairs Q — TL and the nature of the extension TL — > Q for the nine single structures 
investigated in this paper p7| , p6[ |. Out of the seven balanced cases, five have identical point groups and two have 
identical Bravais lattices with P a being the inversion. In the next section, we will discuss our numerical results for 
Fourier series of functions $(r), whose $ = iso-surfaces approximate triply periodic minimal surfaces. However, a 
very simple approximation can be obtained by only considering the space group information. The use of a single Fourier 
star leads to the so-called nodal approximation, which does not need any Fourier amplitude. For balanced structures, 
it is essential to consider the correct black-and-white symmetry. If one Fourier star is not sufficient to represent the 
topology correctly, one has to consider another one and adjust its Fourier amplitude by visual inspection. In the last 
column of Tab. pi, nodal approximations are given for the nine single structures investigated. They have been derived 
in the context of triply periodic zero potential surfaces of ionic crystals [po[ and are widely used as approximations 
for TPMS ( e.g. in Refs. jl2],[36]] ) . However, it is well known that their properties can differ considerably from those of 
real TPMS |37| . In the next section we will show that the quality of the nodal approximation varies considerably for 
different structures. 
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IV. TRIPLY PERIODIC STRUCTURES IN GINZBURG-LANDAU MODELS 



A. Ginzburg-Landau model 

The simplest Ginzburg-Landau model Jl6| , which successfully describes the phase behavior and mesoscopic structure 
of ternary amphiphilic systems, contains a single, scalar order parameter field <fr(r); thus, amphiphilic degrees of 
freedom are integrated out. The model is defined by the free-energy functional 

F[$]= J dr {(A$) 2 + 5 (<i>)(V$) 2 +/(*)} . (7) 

For ternary amphiphilic systems at the phase inversion temperature, the free-energy functional has to be symmetrical 
in water and oil, i.e. invariant under the transformation $ — > — $. A convenient form for the free-energy density and 
the structural parameter of the homogeneous phases are pq | 

/($) = ($ + l) 2 ($-l) 2 ($ 2 + / ) , .g($)=.g +.g 2 $ 2 . (8) 

Here, the three minima at —1, and 1 correspond to excess water, microemulsion and excess oil phases, respectively, 
and fo acts as a chemical potential for amphiphilcs. The behavior of the correlation function ($(ri)$(r2)) is controlled 
by the parameters go and g 2 in the microemulsion and excess phases, respectively; they are therefore related to the 
amphiphilic strength and the solubility of the amphiphile Strongly structured microemulsions are obtained for 
go < 0, while unstructured excess phases require (g 2 + go) to be positive and sufficiently large. 

The elastic properties of amphiphilic monolayers are described by the Canham-H elf rich Hamiltonian |39 40 



H = j ' dA {a + 2k (H - c ) 2 + H K) , (9) 

where the integration extends over the surface, H and K are mean and Gaussian curvature, respectively, and the 
elastic moduli are the surface tension a, the spontaneous curvature Co, the bending rigidity k and the saddle-splay 
modulus R. If the $ — iso-surfaces are identified with the location of the amphiphilic monolayers, the free energy 
of the Ginzburg-Landau model (ffi) can be well approximated by a sum of the curvature energy — calculated from 
Eq. ([)]) — and a direct interaction between the interfaces, which falls off exponentially with distance |l^,^8|. Due 
to the water-oil symmetry the spontaneous curvature Co vanishes. The values of the other elastic moduli can be 
calculated numerically from the profile <f> s (z), which minimizes the free energy of a planar water-oil interface [ jl7| . 

With /o, go and g 2 , the model has a three-dimensional parameter space. We investigate bicontinuous cubic phases 
in the mean-field approximation. For fo < the excess phases W/O are stable, for fo > the microemulsion ME. 
Moreover, for sufficiently negative go and not too large g 2 , a (singly periodic) lamellar phase L a is found to be stable. 
Fig. U shows cuts through the parameter space for go — —3.0 and go — —4.5. At each point of parameter space, the 
elastic moduli can be determined numerically. In the region where a < 0, the system gains free energy by forming 
interfaces; the lamellar phase is stabilized by the short-ranged, repulsive part of the interaction between them. In the 
vicinity of the a — 0-line with a > 0, the lamellar phase still exists — stabilized by the long-ranged, attractive part of 
the interaction between interfaces. However, for a given go there is exactly one point where the phase boundary and 
the a = 0-line touch each other [|dj. Here the interaction between lamellae vanishes and their distance diverges. We 
call this point unbinding point since this divergence resembles a continuous unbinding transition 1 42 4^] . We will later 
make use of the fact that in the vicinity of the unbinding point the free energy of the Ginzburg-Landau model (]?]) 
is dominated by its interfacial contributions. To our knowledge, no other phases are stable in the Ginzburg-Landau 
model. The (doubly periodic) hexagonal phase as well as the (triply periodic) cubic phases investigated below are 
only metastable. If one raises go to go — —2, the lamellar channel vanishes, and microemulsion and excess phases 
coexist for fo = and large g 2 . As go approaches zero from below (with g 2 > 0), the lamellar phase disappears. 



B. Generating triply periodic structures 



In order to generate the structures described in Sec. || as local minima of the Ginzburg-Landau functional, we 
construct the Fourier series for the corresponding space group pair (listed in Tab. ||), as described in Sec. HI, and 
terminate it after N stars. For a given set of parameters (/o, go, 52), the free-energy functional (0) becomes a function 
of the Fourier amplitudes Ai and Bi (1 < i < N), the constant mode Ao and the lattice constant a. Although the 
integral in Eq. (Q) could be solved exactly by using the orthogonality of trigonometric functions, we use Gaussian 



5 



integration which is equivalent in efficiency and accuracy, but easier to implement. For minimization we use conjugate 
gradients and as initial profile the nodal approximations given in Tab. 0. Note that certain structures like P and C(P) 
need the same Fourier ansatz but different initial profiles. In order to obtain initial profiles for the double structures, 
we transform the nodal approximations of the corresponding single structures according to $ — > 2$ 2 — 1 in real space 
and then transform these profiles to Fourier space by using Eq. (|^). In order to test for convergence of the Fourier 
series, we increase N successively and require both the values for the free-energy density / and the lattice constant a 
to level off. With a typical workstation N 100 is feasible, but N sa 30 turns out to be sufficient in most cases. If a 
profile has converged, the Fourier amplitudes fall off exponentially for large |K,-|. 

We have generated the nine single structures discussed in Sec. || and the double structures corresponding to the 
four most relevant single structures. In Fig. || we visualize single and double structures of P, D, G and I-WP by 
drawing their $ = iso-surfaces in a conventional unit cell. The two interfaces of the double structures form a bilayer 
wrapped onto the interfaces of the single structures. For /o = 0.0, go = —3.0 and gi — 7.01 our numerical results are 



summarized in Tab. III. This point is chosen here in order to compare our results with those of Gozdz and Holyst 
p5| . For each structure, we give the free energy density /, the lattice constant a and the volume fraction of oil 
v = 1/V JcZr($(r) + l)/2. The structures are ordered with respect to their free-energy density /; thus, for the single 
structures we find the hierarchy G - S - D - I-WP - P etc. The scaled surface area A/ a 2 and the Euler characteristic \ 
are calculated for the conventional unit cell by triangulating the $ = iso-surface with the marching cube algorithm. 
By successively refining the discretization, we calculate the sum Am over the surface areas of the triangles as a function 
of their number M, and then extrapolate to A — A^ by fitting Am to a linear 1/M -dependence. The values obtained 
for P, D, G and I-WP turn out to be very close to the ones for the corresponding minimal surfaces as calculated from 
the Weierstrass representations and given in Tab. |. The Euler characteristic \ can be calculated as \ = C — E/A 
where C is the number of cubes and E is the number of cube edges cut by the surface. All our numerical results agree 
very well with those reported by Gozdz and Holyst in Ref. [Gl| for G, D, P and I-WP. 

In order to evaluate the curvature properties of the $ = iso-surfaces, we make use of the fact that — due to the 
Fourier ansatz — the field 4>(r) is known analytically. The mean and Gaussian curvatures, @,|37j 



1 1 
H = — V • n , K =- 

2 ' 2 



(V • nr ^ iln, d, 



jit. 



(10) 



where n = V<&/|V$| is the surface normal vector, can then be calculated exactly. In Fig. [| we plot the distribution 
of H and K over the <& = iso-surface as histograms for the structures P, D, G and I-WP in a conventional unit cell. 
Numerical inaccuracies arise here only from the calculation of the position vectors of the $ = iso-surface with the 
marching cube algorithm, and from the triangle areas AAi, which enter as weighting factors of H and K for each 
plaquette of the triangulation. The mean curvature is close to zero everywhere. For the balanced structures P, D and 
G it is distributed symmetrically around H = 0. We conclude that the $ = iso-surfaces of our numerical solutions 
are very close, but not identical to the corresponding TPMS. For the Gaussi an c urvature, we also plot the distribution 
obtained from the Weierstrass representation by numerically evaluating Eq. ( |A3j ). This is done by evaluating K(x+iy) 
on a square lattice in x and y, which covers the fundamental domains for P, D, G and I-WP as given in Appendix [A|, 
and collect the values in a histogram with weights dA{x + iy). The resulting distributions of K are scaled to unit 
lattice constant by using the Weierstrass lattice constants given in Tab. [jj and normalized to J dKp(K) = 1. For P, D 
and G, these distributions p(K) are identical (apart from the scale) since they are related to each other by a Bonnet 
transformation. Fig. |] shows that the numerical distributions agree very well with the exact Weierstrass results. The 
surface area is not very sensitive to the detailed shape of the interfaces in these cases; the values for the full solutions 
and the nodal approximations differ only by a few percent. Finally, we can test the procedure used here by employing 
the Gauss-Bonnet theorem to calculate the Euler characteristic as x = 1/(2tt) / dAK. We find an agreement to three 
relevant digits. 

The ^-distributions of the single structures investigated, for which no Weierstrass representations are known, are 
shown in Fig. ^. These surfaces typically have a more complicated structure within the unit cell, which is reflected in 
both a larger number of peaks and in a larger extremal value of the Gaussian curvature, K m i n . 

The Fourier ansatz yields representations of TPMS which are easy to document and thus straightforward to use in 



further investigations. In Tab. IV, we give improved nodal approximations which consist of up to six Fourier stars. 
When the number of Fourier stars N is increased, the free-energy density / decreases monotonically, but the quality 
of the approximation for the corresponding TPMS does not necessarily improve in the same way. Therefore for each 
structure we choose an optimal value of N . The quality of the approximation is judged from the distribution of H and 
K over the surface, as described above. In Tab. |v|, we compare the nodal approximations from Tab. ||, the improved 



nodal approximations from Tab. IV and the exact Weierstrass representations by monitoring H maxi the maximum of 
\H I on the iso-sur face, y/(H T ), the square of the variance of H on the iso-surface, and K m i n , the minimal value of K 
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on the iso-surface. For the exact Weierstrass representation, one has H max = -J (H 2 ) = 0. For the improved nodal 
approximations, the values for H max and yj (H 2 ) improve by nearly one order of magnitude when compared with 
the nodal approximations. Also, their distributions of K are very close to the ones obtained from the Weierstrass 
representations (the only exception is G, whose nodal approximation is already quite good, in particularly for its 
if-distribution) . Thus by adding just a few more modes, we can considerably improve on the widely used nodal 
approximations. Finally, we want to point out that the values of the amplitudes for P, D and I-WP given in Tab. IV 
are of the same magnitude as those calculated from the representations for TPMS obtained in Ref . p5j . The different 
values arise from the different shapes of the order parameter profile through the interface. The Ginzburg-Landau 
model gives interfaces with a finite width, while the calculation based on TPMS assumes sharp interfaces. In fact, 
the amplitudes of the Ginzburg-Landau model should approach the sharp-interface results as the unbinding point is 
approached. As the interface effectively sharpens, the Fourier amplitudes decay more and more slowly as a function of 
the wave number |K|, with a |K| -2 behavior in the limit of step-like interfaces — in agreement with the well-known 
|K| -4 Porod-law for the scattering intensity of sharp interfaces fiqj . Note that Tab. [V now provides data for G for 
the first time, which in fact is the cubic bicontinuous phase most relevant for amphiphilic systems j(] 



C. Hierarchy of structures 

In order to investigate the relative stability of the cubic structures in the symmetric Ginzburg-Landau model, we 
now consider their interfacial properties. For a triply-periodic cubic structure, the free-energy density within the 
curvature model @ reads 

fcurv = \ M*) + ^{ 2k J dA H2 + 2R7T *) ( U ) 

where A* — A/ a 2 is the scaled area per unit cell; the surface area A, the integration of H 2 and the Euler characteristic 
X once again refer to the conventional unit cell. Both terms in brackets are scale invariant, i.e. they do not depend 
on the lattice constant a. If the elastic moduli are calculated for the points in the phase diagram where we minimize 
the functional, we find a < 0, k > and R < 0. In particular, for /o = 0.0, go — —3.0 and 52 = 7.01 we obtain 
(7 = —0.84587, k = 2.36197 and R = —0.97646. This explains why the bicontinuous phases cannot be stable in 
our model. The negative surface tension favors modulated phases, and the bending rigidity favors minimal surfaces; 
however, only the lamellar phase is not disfavored by the negative saddle-splay modulus. We now also understand 
why the cubic structures are stable in the Ginzburg-Landau model with respect to a variation of a. There is a balance 
between the negative surface tension term, which favors small values of a, and the positive curvature contributions 
which favor large a. The minimization of Eq. (|ll]) with respect to a yields 

_ ( 6TrR X + 6 K fdAH 2 y ( A \ h ( (HO 3 \ k no \ 

amm ~{ \<j\ A * J ' /mm ~ {27) \2Rir X + 2KfdAH 2 ) ' 1 ' 

where A* is assumed to be (approximately) independent of a. We can now understand why the single structures 
are found to be so close to minimal surfaces; the minimum of J dAH 2 (the so-called Willmore problem fl47j]) in this 
case is H = 0, which in turn minimizes f m in- Note that this reasoning is not rigorous, since the minimization of 
the free-energy functional with respect to lattice constant and shape are not independent; the latter step determines 
A/a 2 , which is taken to be constant in the first step. 

We first discuss the single structures without the "complicated" phases C(P), C(D), S and C(Y) — which have 
larger lattice constants, a more pronounced modulation within the unit cell, and stronger interactions between the 
surfaces. Then we numerically find / dAH 2 10 -3 (compare Fig. ||), so that Eq. ( p^ ) becomes 

amm = (rf) ' fmin = " (27) (jk) r (13) 

where T = (A* 3 /2n\x\)i = (A 3 /2ir\ X \a 6 ) 2 is t he t opology index. Its exact and numerical values for the single 



structures investigated is given in Tab. Q and Tab. |III|, respectively. This geometrical quantity is independent both of 
lattice constant and choice of unit cell. It can be considered to be a measure for the porosity of the structure (the 
larger its value, the less porous) as well as for the specific surface area (the larger its value, the more surface area per 
volume). Its relevance for amphiphilic systems is well known fl48| [50|. In fact it follows from the isoperimetric relations 
in three-dimensional space that A 3 /|x|a 6 is the only invariant combination of Minkowski functionals for vanishing 
integral mean curvature KM. 
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A comparison of Eq. (|T^) with our numerical results from Tab. Ill shows that these formulae systematically predict 
a lattice constant and a free-energy density, which are too large and too low by about 20%, respectively. However, the 
hierarchy of structures as predicted by f m % n turns out to be exactly the same as given in Tab. [II by the full numerical 
results for /: G - D - I-WP - P - F-RD. We therefore conclude that the gyroid structure is the most stable structure since 
it has the smallest porosity (the largest topology index). This corresponds to the fact that topologically the gyroid's 
labyrinths have the smallest connectivity of all structures considered — G is the only structure with only three lines 
meeting at one vertex (D and P have four and six, respectively) . With the topology index, we have found a universal 
geometrical criterion for the relative stability of the various single bicontinuous cubic phases, which also resolves 
the debate about the degeneracy of minimal surfaces in the Canham-Helfrich approach for vanishing spontaneous 
curvature. In order to explain the non-degeneracy observed experimentally for binary lipid-water systems, higher 
order terms |32| and frustration of chain stretching ||^3| have been considered. In our description both effects are not 
necessary. In fact, in a (balanced) ternary system, chain stretching does not provide a plausible mechanism for lifting 
the degeneracy of the free energies of different TPMS, since the presence of oil relieves the frustration in the chain 
conformations. In the part of the phase diagram, where ordered phases are stable, the free-energy density should 
contain a negative surface tension contribution (compare Refs. |^6| , p8| ) and the topological term — which due to the 
Gauss-Bonnet theorem is often neglected, but in fact prevents the lattice constant from shrinking to zero Q|. The 
same reasoning might be applied to binary systems by identifying monolayers with bilayers and oil and water with 
water I and water II, respectively. Although the presence of a negative surface tension is essential in our argument, 
we want to emphasize that its magnitude can still be very small. 

There are two main reasons why Eq. (|l^) yields the correct hierarchy in / m j n , but does not reproduce the numerical 
values very well. First, by focussing on the interfacial properties, we have neglected the contributions to the free- 
energy density due to direct interactions between the interfaces, and second, the reasoning leading to Eq. (|l2|) is not 
rigorous. In order to check whether this is indeed the origin of the numerical discrepancy, we can make use of the 
existence of the unbinding point in the phase diagram, Fig. |^. At this point, both a and the interaction vanish and 
the structures investigated should become exactly minimal surfaces. However, at the same point the lattice constant 
diverges, and both the minimization and the Fourier ansatz become unfeasible. Therefore, we study the approach to 
the unbinding point along the path shown in Fig. ^. Numerically we find for P, D and G that the deviation of the 
values for / (a) from Eq. (12) relative to our numerical results reduces from 14% (20%) through 11% (16%) to 5% 
(7%) for the three sets of parameters considered. This demonstrates the convergence towards minimal surfaces, and 
verifies our reasoning above. 

In contrast to the case of "simple" single structures discussed so far, the topology index T does not predict the 



hierarchy found numerically for the more complicated single structures C(P), C(D), S and C(Y), compare Tab. Ill 



This may be due to numerical uncertainties in the value of T, since no exact results are available for the surface area 
A/ a 2 . However, both the low Ginzburg-Landau free-energy density and the large topology index of the S-surface 
indicate that this structure should be comparable in stability with the G-surface. We believe that double structures 
are not favored in our model because their interfaces are not minimal surfaces, and can be better described as surfaces 
of constant mean curvature. 

One of the intriguing aspects of the Ginzburg-Landau model is that it has a very rugged energy landscape with many 
local minima. In fact, we were able to find many more interesting structures, including more complicated minimal 
surfaces and surfaces which contain both saddle-shaped pieces and pieces with positive Gaussian curvature. The latter 
result is related to the observation that we essentially solve the Willmore problem, which allows for solutions which 
have both regions with K > and regions with K < — for example, the Clifford torus p7|. 



V. SCATTERING INTENSITIES AND NMR-SPECTRA 



We expect the nine single and four double structures investigated in this paper to include all physically relevant 
bicontinuous cubic phases in ternary systems near their phase-inversion temperature. The representations of these 
phases can now be used to simulate data for those experimental techniques, which mainly depend on the geometry of 
the structures. This includes small angle scattering (SAS), transmission electron microscopy (TEM) and nuclear mag- 
netic resonance (NMR). In particular, our representations can be used to analyze experimental data for mesoporous 
systems when bicontinuous cubic phases have been used as templates [ [l3| . 

Since their typical length scales are in the nanometer range, amphiphilic structures can be investigated by X-ray 
and neutron small angle scattering (SAXS and SANS). The scattering function S(k) is obtained from the Fourier 
transform $(k) of the density <E>(r) of scattering lengths as S(k) = |$(k)| 2 . For SANS, the scattering contrast can 
be varied by deuterating the sample. For balanced single structures, it is clear from above that with bulk contrast 
one measures space group 7i, with film contrast space group Q . Thus, by varying the contrast in SANS-experiments, 
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the space group pair of the investigated structure can be determined |3^|35|] . With our representations of the cubic 
structures, it is straightforward to calculate the corresponding scattering amplitudes which are needed to analyze the 
experimental data. 

It is quite clear from our results why the program to identify a structure from its scattering intensity alone is 
rather difficult. For the single structures, the Fourier amplitudes decay so rapidly with increasing wave vector that 
the scattering intensity is dominated by the first peak, while the higher peaks are hardly detectable. Even for double 
structures (or film contrast for single structures in SANS) the situation hardly improves, since the scattering intensity 
also does not feature more than two or three relevant peaks. Note that even if the correct space group were extracted, 
one still could not be certain about the type of minimal surface; G and S, for example, have the same space group 
laid (No. 230) when measured in film contrast. 

It has been pointed out by Anderson [^8[ that deuterium NMR is another experimental technique which yields char- 
acteristic fingerprints of bicontinuous cubic phases. Three conditions have to be met in this case; first, the amphiphilc 
has to be deuterated and distributed uniformly over the surface, second, the sample has to be monocrystalline, and 
finally, the surfactant sheets have to be polymerized in order to suppress diffusion within the surface. Since deuterium 
has a nucleus with 7 = 1 and a small magnetic moment, the quadrupolar interactions dominate and the dipolar ones 
can be neglected. The 2 if -NMR experiments then essentially measure the distribution f(x) of x = (3cos 2 p — l)/2, 
with — 4 < x < 1, where p is the angle between the external magnetic field and the amphiphilic director, which cor- 
responds to the normal vector of the interface [57[ ]. In fact, the electrostatic field at the nucleus does not distinguish 
between x and —x and the bandshape is g(x) — f(x) + f(—x), with —1 < x < 1. 

For magnetic fields parallel to the (lOO)-direction, f(x) can be calculated exactly for P, D, G and I-WP from their 
Weierstrass representations. The calculation for P, D and G has been performed by Anderson |5q| . We introduce 
spherical coordinates (0, ip) with the polar axis in the direction of the magnetic field, so that p — 9. The Weierstrass 
representation uses the complex w-plane (see Appendix ^), which is parametrized in polar coordinates (r, i^). At every 



point bj, the normal vector follows by stereographic projection to the unit sphere. Eq. (A2) then gives the one-to-one 
correspondence 

r 4 - 4r 2 + 1 , , 

X = — ^r— . (14) 

(1 + r 2 ) 2 v ; 



With r' = r(x') and Eq. (A3), the distribution becomes 
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dy\R{r'^)\ 2 (15) 
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where R is the generating function given for P, D, G and I-WP in Appendix We denote the remaining integral by 
I(r'). P, D and G yield the same result [ff], 

32 K(k) 

IpDG{r) = \ = , (16) 

r 4 V(ei - l)(e 2 + 1) 



since they are related by a Bonnet transformation. Here, K{k) is the complete elliptical integral of the first kind, and 

2( £l - eg) 
( £l - l)(e 2 + 1) 



a+ + a_r 8 a_ + a+r 8 7 ± ^48 ;2 2(ei - e 2 ) 

e i = 1 ' £ 2 = 2 ' a ± = — n — i k = 7 T77 — rrr • ( 17 ) 



For I-WP, we find from Eq. (0) 



( r 2 _ r iO)f s \ l-r 



where P v {x) is a Legendre function of the first kind. Two values of r correspond to any given x in Eq. (|lj), < r < 1 
for the southern hemisphere and 1 < r < oo for the northern hemisphere. However, for the direction and the structures 
considered both x and dA do not change when the northern is mapped onto the southern hemisphere by a reflection 







through the xy-plane (this can be verified by inverting the complex plane in the unit circle). It is therefore sufficient 
to evaluate Eq. (|l5| ) for < r < 1. Then the bandshapes g(x) = f(x) + f(—x) for both P, D, G and I-WP are found 
to have peaks at x = ±i, which are already present in the Pake pattern, f(x) = 2n/ \/6x + 3, of isotropic samples. 
However, additional peaks appear at x = for P, D and G and at x — ±1 for I-WP, which correspond to the positions 
of the flat points ((111) for P, D and G, and (100) for I-WP). 

For the same structures, but other directions, and for other structures, an exact calculation is not possible. However, 
we can numerically determine the distribution f(x) for any direction and any structure in the same way used above 
to calculate the distribution of H and K over the surface. In fact, with the Fourier ansatz, this can be done with 
much better resolution than for a real-space representation [^9) . In Fig. ^ we show our results for P and I-WP and 
the three high-symmetry directions. For the (lOO)-direction, we also show the exact results following from Eq. jl5| ) 
with Eq. ([16]) and Eq. ([18]) for P and I-WP, respectively; the agreement with the numerical results is excellent. In 
Fig. [t] our results are shown for C(P) and F-RD for which no Weierstrass representations are known. 

VI. SUMMARY 

We have presented here a systematic investigation of bicontinuous cubic phases in ternary amphiphilic systems. 
In these structures the amphiphilic monolayers form triply periodic surfaces with cubic symmetry. We distinguish 
between single and double structures which are characterized by one or two monolayers, respectively. In order to 
further classify the structures within each of these groups, we used the crystallographic classification of triply periodic 
minimal surfaces by Fischer and Koch. Thus for single structures we distinguish between balanced and non-balanced 
structures; the first class is further subdivided into spanning and non-spanning structures. Finally we end up with a 
list of nine single structures of interest (P, C(P), D, C(D), S, C(Y), G, I-WP, F-RD). For each of these there exists a 
corresponding double structure from which we considered those corresponding to the simple and physically relevant 
single structures P, D, G and I-WP. 

In the framework of a Ginzburg-Landau model for ternary amphiphilic systems, we generated the nine single and 
four double structures using the Fourier ansatz and the theories of space groups and color symmetries. Compared to 
real-space minimization, the Fourier approach has the advantage of efficient numerics and easy documentation. For 
P, D, G and I-WP, we gave improved nodal approximations which give much better approximations to triply periodic 
minimal surfaces than the widely used nodal approximations by von Schnering and Nesper. We showed that the free- 
energy density of the single structures can be calculated from an effective surface Hamiltonian with negative surface 
tension, positive bending rigidity and negative saddle-splay modulus. Due to the water-oil symmetry of the model, the 
spontaneous curvature vanishes and structures with zero mean curvature are favored. A comparison with the exact 
Weierstrass representations for P, D, G and I-WP shows that the single structures can be made to closely approach 
triply periodic minimal surfaces by appropriately tuning the model parameters. For vanishing mean curvature term, 
their relative stability is determined by the topology index T = (A 3 /27r|x|a 6 ) 5 . This explains the hierarchy G - D - 
I-WP - P. Thus for water-oil symmetry the single gyroid is the most stable cubic bicontinuous phase since it has the 
smallest porosity. 

The representations obtained for both single and double structures can now be used for further physical investiga- 
tions. We have employed them to calculate the distribution of the Gaussian curvature for C(P), C(D), S, C(Y) and 
F-RD, which were not known before. Furthermore, we have determined several quantities which can be measured in 
SAS- and 2 i?-NMR-experiments. 

In our Ginzburg-Landau model for balanced ternary systems, modulated phases are favored whose interfaces are 
minimal surfaces. However, only the lamellar phase is stable, since the bicontinuous phases are disfavored by the 
negative saddle-splay modulus. Thus, other energetic contributions have to be considered to stabilize bicontinuous 
phases. This can be long-ranged interactions (like the van-der-Waals or electrostatic interactions) or higher order 
terms in the curvature energy. Assuming that these contributions typically have a similar effect on all bicontinuous 
phases, we then expect the gyroid phase G to be most prominent. One of the intriguing conclusion of our calculations 
is the large relative stability of the S-surface, as indicated by its low free-energy density in the Ginzburg-Landau 
model and by its large topology index. Thus, for single cubic phases in balanced ternary amphiphilic systems, the S- 
surface could be a possible alternative to the G-surface, whose double version seems to be so ubiquitous in lipid-water 
mixtures. 

We want to emphasize that bicontinuous phases can of course be induced by a positive saddle-splay modulus. 
However, in this case the lattice constant should be on the order of the size of the amphiphilic molecules, and the 
application of the curvature energy model becomes questionable. This may indeed be the case in many lipid-water 
systems, where the amphiphile volume fraction in the bicontinuous cubic phases is larger than 50%. 
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APPENDIX A: WEIERSTRASS REPRESENTATIONS 




Exact representations are known for four cubic TPMS |2l],^2|,^0|,^l| . Consider the composite mapping of the surface 
into the complex plane, which consists of two parts; first, each point of the surface is mapped to a point on unit 
sphere which is defined by the normal vector on the surface, then the unit sphere is mapped into the complex plane 
by stereographic projection. The Weierstrass representation formulae invert this mapping and therefore map certain 
complicated regions ft C C of the complex plane onto a fundamental piece of the TPMS, 

i(x,y)=Rej^ dz R(z) + z l ) \ . (Al) 

where to = x + i y. The replication of this surface segment with the symmetries of the space group TL of the 
oriented surface then gives the whole TPMS. It can be shown that each meromorphic function R(uj) corresponds to 
a minimal surface. However, only few meromorphic functions lead to embedded TPMS. For the D-surface, one has 
R(u>) = (lu 8 — 14a; 4 + 1) _ 2 ; and f2 is the common area of the four circles with radii V% around the points (±1 ± i)/ \[2. 
The Bonnet transformation R{io) — > e ie R(u>) with 9 = it/2 transforms the D-surface piece into one for the P-surface, 
and 9 = 38.015° does the same for the G-surface. The generating function for I-WP was found only recently to be 
R(lu) — (uj(lu 4: + |6^j6^| . Here fl is the common area of the circle with radius \f2 around the point (1 + i)/V2 

with the unit circle, the lower half plane and the lower half plane rotated anti-clockwise by ir/4. 

Given a generating function R(oj) for a certain TPMS, several of its properties can be readily calculated. At a given 
point u> = x + iy = rexp(i(^), the normal are determined by stereographic projection from the complex plane to the 
unit sphere, 

cos9=r ^j\ ( A2 ) 
in spherical coordinates {9,ip). Gaussian curvature and differential area clement follow as 

K ^ = I p/ , 2 ," 4 m , 2vt . dA{u) = \R{uj)\ 2 (1 + M 2 ) 2 dxdy . (A3) 

Therefore the poles of R{ui) correspond to the (isolated) flat points of the surface (that is points with K = 0). Since K 
and dA depend only on the modulus of R, surfaces related by a Bonnet transformation have the same distribution of 
K and the same surface area for one fundamental domain. Tab. || shows data which are obtained from the Weierstrass 
representations . 
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N d 


X 


a 


A* 


r 


V 


p 


12 


-4 


2.1565 


2.345103 


0.716346 


0.5 


D 


48 


-16 


3.3715 


3.837785 


0.749844 


0.5 


G 


24 


-8 


2.6562 


3.091444 


0.766668 


0.5 


I-WP 


48 


-12 


7.9499 


3.464102 


0.742515 


0.536 



TABLE I. Properties of the four cubic triply periodic minimal surfaces G, D, I-WP and P, for which (exact) Weierstrass 
representations are known. Nd is the number of copies of the fundamental domain which are necessary to build up the surface 
in one conventional unit cell. \ is the Euler characteristic, a the lattice constant, and A* = A/a 2 the scaled surface area in the 

., r. 1 

conventional unit cell. V — (A /2ir\x\a )' is the topology index and v the volume fraction of one of the two labyrinths. 





n 


g 


n^g 


nodal approximations 


p 


Pm3m (221) 


ImZm (229) 


p -> i 


Eccc(100) 


C(P) 


Pm3m (221) 


Im3m (229) 


p -> / 


Eccc(100) + Eccc(lll) 


D 


Fd3m (227) 


PnZm (224) 


F ^ P 


(Eccc + Ecss + Escs + Esse) (111) 


C(D) 


Fd3m (227) 


Pn3m (224) 


F -> P 


(Eccc + Ecss - Escs - Essc)(311) 


S 


/43d (220) 


Ia3d (230) 


43m — ► m3m 


(Ecsc + Occs)(211) 


C(Y) 


P4 3 32 (212) 


J4i32 (214) 




- (Eccc + Osss + Esss + Occc)(lll) 
+ 6 (Esse + Oscc + Eccs + Ocss)(210) 


G 


J4i32 (214) 


IaM (230) 


432 — > m3m 


(Escc + Ocsc)(110) 


I-WP 


Im3m (229) 






2 Eccc(llO) - Eccc(200) 


F-RD 


Fm3m (225) 






4 Eccc(lll) - 3 Eccc(220) 



TABLE II. Group - subgroup pairs g — Ti, the relation 7i — » g between them and nodal approximations for the cubic 
minimal surfaces. They are grouped into three classes: spanning with three-dimensional nets, balanced but non-spanning, 
and non-balanced, g and H for balanced single structures differ either in Bravais lattice or in point group (but not in 
both). Non-balanced single structures have TL = g. The Fourier ansatz has to be chosen as follows: TL with color sym- 
metry (given by H — > g) for single structure, and g for double structure. Nodal approximations only consider the space 
group information given by g — Ti. We use the definitions Epqr = p(hx)q(ky)r(lz) + p(hy)q(kz)r(lx) + p(hz)q(kx)r(ly) and 
Opqr — p(hx)q(kz)r(ly) + p(hy)q(kx)r(lz) + p(hz)q(ky)r(lx), where p, q and r can be either c (cosine) or s (sine). 
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f 
J 






A* 


v 

A 




G 


-0 19096 


10.08 


0.500 


3 09140 


-8 


76665 


s 


-0 18962 


17.72 


0.500 


5.41454 


-40 


0.79474 




-fl i 8870 


12.56 


0.500 


3 83755 


-16 


0.74978 


I-WP 


-0.18112 


11.82 


0.527 


3.46367 


-12 


0.74238 


P 


-0.18109 


7.89 


0.500 


2.34516 


-4 


0.71637 


C(Y) 


-0.18061 


14.92 


0.500 


4.46108 


-24 


0.76730 


C(D) 


-0.17382 


28.93 


0.500 


8.25578 


-144 


0.78862 


F-RD 


-0.16311 


17.36 


0.532 


4.75564 


-40 


0.65417 


C(P) 


-0.16239 


14.13 


0.500 


3.80938 


-16 


0.74154 


GG 


-0.18850 


18.47 


0.495 


5.33712 


-16 




DD 


-0.18549 


11.59 


0.502 


3.30239 


-4 




PP 


-0.17700 


14.98 


0.523 


4.09739 


-8 




IWP2 


-0.16193 


21.55 


0.532 


6.29274 


-24 





TABLE III. Properties of cubic bicontinuous phases as calculated numerically from the Ginzburg-Landau functional at the 
point /o = 0.0, go = —3.0 and g2 = 7.01. The lamellar phase is stable at this point with / = —0.20807, the hexagonal phase 
metastable with / = —0.17552. The nine single (top) and four double (bottom) cubic structures are ordered with respect to 
their free-energy density /. a is the lattice constant, v the volume fraction of oil. The values for scaled surface area A* = A/a 2 
and Euler characteristic \ are given for the conventional unit cell. From these numbers the topology index Y can be calculated. 
Note that the symmetry between water and oil is broken not only for non-balanced, but also for double structures. At this 
point of the phase diagram, about 50 Fourier modes are sufficient and have been used for all structures. 





mode 


function 


amplitude 


p 


(1 0) 


Eccc -(- Occc 


0.2260 




(1 1 1) 


Eccc + Occc 


-0.0516 




(2 1 0) 


Eccc + Occc 


-0.0196 




(3 0) 


Eccc + Occc 


-0.0027 


D 


(1 1 1) 


Eccc + Ecss + Escs + Esse + Occc + Ocss + Oscs + Ossc 


0.1407 




(3 1 1) 


Eccc + Ecss - Escs - Esse + Occc + Ocss - Oscs - Ossc 


0.0209 




(3 1 3) 


Eccc - Ecss + Escs - Esse + Occc - Ocss + Oscs - Ossc 


-0.0138 




(1 1 5) 


Eccc + Ecss + Escs + Esse + Occc + Ocss + Oscs + Ossc 


-0.0028 




(3 3 3) 


Eccc + Ecss + Escs + Esse + Occc + Ocss + Oscs + Ossc 


-0.0021 




(3 1 5) 


Eccc + Ecss - Escs - Esse + Occc + Ocss - Oscs - Ossc 


0.0011 


G 


(1 1 0) 


Escc + Ocsc 


0.3435 




(0 3 1) 


Ecsc - Occs 


-0.0202 




(2 2 2) 


Esss + Osss 


-0.0106 




(2 1 3) 


Ecsc + Occs 


0.0298 




(3 3) 


Eccs + Oscc 


0.0016 




(1 1 4) 


Escc + Ocsc 


-0.0011 


FWP 


(0 0) 


1 


0.0652 




(1 1 0) 


Eccc 


0.5739 




(2 0) 


Eccc 


-0.1712 




(2 1 1) 


Eccc 


-0.1314 




(2 2 0) 


Eccc 


0.0184 



TABLE IV. Improved nodal approximations for the four most relevant cubic structures P, D, G and I-WP. In contrast to 
the nodal approximations given in Tab.Jll|, the extremal values of the functions 3>(r) are close to ±1, since they are calculated 
from the Ginzburg-Landau model ([?]), Epqr and Opqr, with p,q,r £ {c, s}, are defined in Tab. 



1G 





Nodal 


Improved nodal 


Weierstrass 








K 


Umax 


V( H ' 2 ) 


K 


K 


p 


1.282 


0.579 


-39.48 


0.157 


0.071 


-16.84 


-18.51 


D 


0.384 


0.184 


-39.48 


0.048 


0.022 


-46.22 


-45.24 


G 


0.208 


0.102 


-29.60 


0.078 


0.046 


-29.42 


-28.08 


I-WP 


3.224 


1.326 


-87.34 


0.273 


0.113 


-45.47 


-44.37 



TABLE V. Comparison of nodal approximations, improved nodal approximations and exact Weierstrass representations. 
Hm ax is the maximum of \H\, K m i n the minimum of K, and (H 2 ) = 1/Aj dAH 2 the variance of H over the surface. The form 
of nodal and improved nodal approximations is listed in Tab. H and Tab. |lV|, respectively. 
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D C(D) 

FIG. 1. Different TPMS with the same space group and the same system of straight lines. Spanning minimal surfaces are 
always balanced, thus they are characterized by a group-subgroup pair Q — TL. For P and C(P) this is Im3m (No. 229) - PmZra 
(No. 221) and for D and C(D) this is Pn3m (No. 224) - Fdlm (No. 227). 
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'0 



FIG. 2. Phase diagram for the Ginzburg-Landau model (Q), for (a) go = —3.0 and (b) go = —4.5. Stable phases are the 
microemulsion ME, the lamellar phase L a and the water and oil excess phases W and O. First and second order transitions 
are drawn by solid and dashed lines, respectively. Tricritical points are marked by squares and unbinding points by triangles. 
The dotted line in the L a stability region is the a = 0-line. Metastable cubic phases are calculated at the points marked by 
full circles. 
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FIG. 4. Distribution of mean curvature H (top) and Gaussian curvature K (bottom) in a conventional unit cell, plotted 
as histograms for the structures P, D, G and I-WP, for which Weierstrass representations are known. Radii of curvature are 
measured in units of the lattice constant. For K, we also plot as solid lines the distributions calculated from the Weierstrass 
representations. Due to the existence of Bonnet transformations between P, D and G, they have the same shape but different 
scales in these cases. The parameters in the Ginzburg-Landau model are the same as in Fig. ^a and the /^-distributions are 
normalized to J dKp(K) — 1. For these "simple" structures, about 30 Fourier modes are sufficient. 
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FIG. 5. Distribution of Gaussian curvature K for single structures, for which no Weierstrass representations are known. The 
normalization is the same as in Fig. ^[ These "complicated" structures have more pronounced modulation within the unit cell 
and larger lattice constants, therefore we use about 50 to 70 Fourier modes for the numerical minimization. 
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FIG. 6. 2 i/-bandshapes for the structures P and I-WP and for the three high-symmetry directions (100), (110) and (111). 
Results from the Weierstrass representations for the (lOO)-direction are drawn as solid lines. 
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FIG. 7. 2 i/-bandshapes for the structures C(P) and F-RD and for the three high-symmetry directions (100), (110) and (111). 
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